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Abstract: 

The possibility of formation of quasibound states and ordered structures in dense 
plasma is investigated. The effective potentials of dense plasma are used. On the basis of 
these models the condition of the ordered structures formation in the system is obtained. 
It is shown that in the dense classical plasma the quasibound states form whereas the 
ordered structures form in the dense semiclassical plasma. 
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I. INTRODUCTION 



It is well known 1-4 that a classical Coulomb system form several ordered structures. As 
shown in Refs. 1 and 2, the strongly coupled one-component plasma will crystallize into a 
centered cubic form at the values of the coupling parameter r> 1. Here, T = e 2 / (ak B T); 
e is the electrical charge; a = (3/47rn) 1//3 is the average distance between particles (Wigner 
- Seitz radius); n is the number density of particles; k B and T are the Boltzmann constant 
and the temperature of plasma, respectively. 

In previous papers 3 ' 4 , we studied the structural and thermodynamic properties of the 
two-component model plasma. The formation of the "near ordering" and anomalous 
increasing of the particle correlation radius as a consequence of microstructure changes 
in plasma have been demonstrated by Monte Carlo simulation method. The "hexatic"- 
like structures which characterized of liquid crystals were found in Refs. 5,6. In Refs. 
7,8 the properties of the shell structures of finite one-component plasma clouds have been 
investigated. The study of ordered structures is especially important in connection with 
the rapid development of dusty plasma. The ordered structures that form in a strongly 
coupled plasma have been first analyzed theoretically in Ref. 9 and have come to be 
known as plasma crystals. The plasma crystals of dust particles were found also in a 
partially ionized plasma in a radio- frequency discharge 10 ' 11 , in the DC glow discharge 12 
and in the UV - induced dusty plasmas under microgravity 13 . 

It is now recognized that the plasma crystals are new material for future technologies. 
Therefore, the investigation of several ordered structures in Coulomb systems plays an 
important role. 

II. INTERACTION MODELS 

In this work we consider a fully ionized, dense (classical and semiclassical) hydrogen 
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plasma. The number density is considered in the range n = n e = rii ~ (10 19 2 x 10 25 ) 
cm~ 3 , and the temperature domain is T ~ (5 x 10 4 -r- 10 6 ) K. 

According to Ref. 14, we can separate a two types of dense plasma. For example at 
9 3> 1 we have a classical dense plasma. Here 9 = k B T/E F , k B denotes the Boltzmann 
constant, T, Ep are a temperature and the Fermi energy, respectively. In this case, the 
quantum mechanical diffraction and symmetry effects are negligible except in short-range 
collisions. When 9 < 0, 1, the electrons are in the state of full Fermi degeneracy. In the 
intermediate region between these cases (9 < 1), the electrons are partially degenerate 
and we have a dense semiclassical plasma. 

Due to the well known long range character of the Coulomb interaction between par- 
ticles in plasma, the correlation effects play an important role for dense plasma. Conse- 
quently, for the dense plasma, the simultaneous interaction of a great number of particles 
should be taken into account. In Ref. 15, an integro-differential equation for the effec- 
tive pair potential has been derived on the basis of Bogolyubov's chain equation. This 
equation takes into account simultaneous correlations of N particles. In the case of three 
particle approximations the expression for the effective potential for dense plasma is 15,16 : 

1 } R l + c( 7 ) ' {) 

Here f(R) = (e~^ R — 1)(1 — e~ 2R )/5 and R = r/r D , where r D is the Debye screening 

length. The potential is expressed in terms of the thermal energy, ^{R) = ^(R)/ksT, 

and 7 = e 2 / {rr,k,BT) is a nonideality plasma parameter. 0(7) is the correction coefficient 

for differrent values of 7: c( 7 ) = -0.008617 + 0.455861 7 - 0.108389 7 2 + 0.009377 7 3 . 

For the semiclassical plasma, the effective potential of Kelbg-Deutsch-Yukhnovskii 17 ' 18 
is usally used: 



*a9(r) = [1 - exp(-r/A Q/3 )] + (2) 
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S a /3S ea k B Tlia(2) exp 



7rln(2)A2 

where \ a p = /i/(27r/x Q , / 3fcsT) 1//2 is the thermal de Broglie wavelength; /i Q( g is the reduced 
mass of electrons and ions. Notice that the effective potential (0) does not account for 
screening effects in plasma and behaves like the Coulomb potential for r — > oo. There- 
fore, in Ref. 19, the effective semiclassical potential was obtained by applying the spline 
approximation to potential (Ef) and the numerical solution of following Equation 19 : 



A# - 3r^ = ±3r^ 2 , (3) 

with the boundary conditions 

* |^o= V/R ; * |^oo= (4) 

at the intersection point. 

The effective potential ^{R) is expressed in units of k%T and R = r/a; A is the 
Laplace operator. In Equation (||) the minus and plus signs correspond to the interaction 
of particles with equal and opposite charges, respectively. The spline approximation was 
performed at the intersection point of Eq.(^) and the numerical solution of © _ ©• This 
potential contains quantum diffraction and symmetry effects at short distances as well as 
screening effects for large distances (see Fig.l). 

III. ORDERED STRUCTURES IN DENSE PLASMA 



Following 20 , let us introduce a potential function for the mathematical description of 
the system state. This function takes into account the collective interaction between the 
components of statistical system: 

U(r,v,t) = [$(r,v,r\v't)f(r\v't)dr'dv'. (5) 

Here $(r, v, r v't) is the potential describing the interaction of particles in the system; 
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a, v are a coordinate and a velocity of particles, respectively. 



If many-particles effects are considered (e.g. three particle effects), the expression for 
U (r, v, t) becomes non-linear: 



U(f,t) = J ^iz(\r — f'\)f(f',v,t)df'dv + J J $123(1^* — r '\ , \f — r "\ , \r ' — r "\) (6) 

f(f v', t)f(r", v", t)df'dv'f"dv". 

In the case of pair central interaction between particles we can written as 

<&(r,v,r',v't) = <&(|r- r'|). (7) 

Taking into consideration Eq.(^), the expression for potential function is obtained: 

U(f,t) = f <S>(\r-r'\)f(r',v,t)dr'dv. (8) 



Considering the effective potentials and the field of collective interactions, we can 
obtain the following set of equations for the potential function U and plasma particles 
density function p 20 : 

U(f) = C J § ef} (\r-r'\)e- u{f " )le dr' (9) 
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p(r) = C exp 
where 9 = 1/ksT] C is the constant. 

We consider 20 the occurrence problem of space-periodical solution for Eq.(^). In this 
case, the expression for U(r) can be written in the following form: 

tp{f) = \J $ e// (|r-r'|)e^'W', (10) 

where tp = —U/9; A = —C/9. Let us that at A = Ao and (p = (po we have a violation of 
space-uniform solution. Let X = X + e, if = (po + u(r) and the equation for u(r) is written 
as: 

^0 + u(r) = (A + e) [ $ e// (|f - f'\)e^ +u ^df', (11) 



After expanding the element of integration as a power series in u(f), we have the following 
expression: 

L[u{f)\ = u{r) -A* J $ e //(|f- r'\)u{f')dr\ (12) 

where X* = X e^°; <x(0) = 4vr / °° ^(s)s 2 ds. 
We consider the following linear equation 

u(r) = A e w J $ e //(|r - f '|)w(f ')dr '. (13) 

According to the occurrence problem of periodical structures from homogeneous medium, 
the solution for u(r) is sought as u(r) = Ae tkr . After inserting this solution into the Eq. 
([13|) and integration, we have the following expression 

e*** = A e^ • • 4tt f°° ^ eff (s)^^-s 2 ds. (14) 

J ruS 

Consequently, in order that ordered structures in the system can form, it is necessary 
that the periodical solution of the following equation of nonlocal statistical mechanics to 
exist 20 : 

e a f°° * / n sin(fcs) 2 , , n , ir . 

- = -4tt / $eff(s) r — ~ s ds = a(k). (15) 

p Jo ks 

In the present work, the behavior of the function a(k) is investigated for the above men- 
tioned models of dense plasma; here k is a wave number. 

The model of dense classical plasma is characterized by the fact that for the low limit 
of integral fllED we must take a certain non-zero value of distance Rq = r/ro- In this case, 
i?o is minimal distance between particles. In Figure 2, the results of numerical analysis 
of expression (|T5| ) on the basis of effective potential (0) are shown. Notice that, with 
increasing of 7, on the dependencies of cr{k) some non-monotonic behavior are observed. 
In this case, the increase of 7 is caused by the density increase (the decreasing of ro). It 
is well known that for a dense plasma the number of particles in Debye sphere Nd < 1. It 
is obtained also that functions cr(fc) form local minimums when Rq ~ 1 (r pa rp). These 
facts can be interpreted as the formation of quasibound states in dense classical plasma. 
The term "quasibound" states 21 ' 22 denotes that the energy e of electron-proton pairs is 
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—e 2 /a < e < 0. Because the size of " quasibound" states is greater than a, it is considered 
that the electron moves over a part of elliptical orbit in the field of the nearest proton, 
thereupon the electron passes into the field of another proton, and so on. 

For the model of dense semiclassical plasma, the problem of introduction of minimal 
distance is absent because the pseudopotential (full line in Fig.l) is limited for r = and 
it is strongly screened when r — > oo. The results of numerical analysis of Eq. are 
represented in Figures 3,4 and 5 for different values of the density parameter r s = a fas 
and the coupling parameter T = Z a Zpe 2 / (aksT). Here, a is the average distance between 
particles and as is the Bohr radius. It is shown that no oscillations on dependencies of 
a ee (k) for electron-electron component are observed. It is probably connected with the 
fact that electrons do not form complexes (structures) between themselves because of 
their little masses. The functions a ep and a pp have pronounced oscillations that can be 
considered as formation of ordered structures. 

In Refs. 2,23 it has been demonstrated that in the intermediate region (1 < F < 150) the 
plasma can be considered as a ordered structures of ions. According to this model, ions 
form a cubic-like space structures due to the strong interactions between particles. Since 
we have a chaotic motion in the system, these structures are realized "on the average". 
Therefore, the energy of this structure is greater then Madelung's energy of simple centered 
cubic lattice. These configurations are called as "ordered structures". 

Consequently in dense semiclassical hydrogen plasma, the "picture" of particles dis- 
tribution can be as follows. Protons are distributed with a certain order and they are 
surrounded by electron clouds. The "proton-electron" system can be considered as a 
quasi-particle. 

This work was supported by Ministry of Sciences of Kazakhstan, Grant MN-2.1/98 
and by INTAS-RFBR, Grant No. 95-1335. Our thanks are due to Yu.L.Klimontovitch, 
V.E.Fortov and A.P.Nefedov for many stimulating discussions. 
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Figure captions 



Figure 1: Effective electron-proton potential for dense semiclassical hydrogen plasma at 
r s — 1 and T = 2 (solid line). Triangles denote the numerical solution of Eq. (^) for a 
dense, classical plasma 15 which accounts for higher-order screening effects. The dashed 
line denotes potential 17 ' 18 which shows quantum corrections at short distances. The dot- 
dashed line: DH potential. 

Figure 2: Reduced a* = a / '(rfjksT) functions for dense classical plasma at Rq = 2. 

Figure 3: Reduced a* e = a ee / (a 3 k B T) functions between electrons of dense semiclassical 
hydrogen plasma at r s = 5. 

Figure 4: Reduced a* ep = a ep / (a 3 kBT) electron-proton functions of dense semiclassical 
hydrogen plasma at r s = 2. 

Figure 5: Reduced a* p = a pp /(a 3 kBT) functions between protons of dense semiclassical 
hydrogen plasma at r s = 1. 



9 




T.S.Ramazanov, M.A.Bekenov, N.F.Baimbetov 
Figure 1: Effective electron-proton potential for dense semiclassical hydrogen plasma at r s = 
1 and r = 2 (solid line). Triangles denote the numerical solution of Eq. (|3|) for a dense, 
classical plasma 15 which accounts for higher-order screening effects. The dashed line denotes 
potential 17,18 which shows quantum corrections at short distances. The dot-dashed line: DH 
potential. 
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T.S.Ramazanov, M.A.Bekenov, N.F.Baimbetov 

Figure 2: Reduced a* = a / '(r^ksT) functions for dense classical plasma at Rq = 2. 



11 




12 




T.S.Ramazanov, M.A.Bekenov, N.F.Baimbetov 

Figure 4: Reduced a* ep = a ep / (a 3 fesT) electron-proton functions of dense semiclassical hydrogen 
plasma at r s = 2. 
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T.S.Ramazanov, M.A.Bekenov, N.F.Baimbetov 

Figure 5: Reduced a* p = a pp /(a 3 kBT) functions between protons of dense semiclassical hydro- 
gen plasma at r s = 1. 
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